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Abstract 

We illustrate an alternative derivation of the viscous regularization of the dif¬ 
fusion equation which was studied in [A. Novick-Cohen and R. L. Pego. Trans. 
Amer. Math. Soc., 324:331-351]. We provide an alternative proof of existence of 
solutions, based on the Galerkin method and on compactness arguments. In ad¬ 
dition, we propose a “non-smooth” variant of the viscous regularization which we 
believe may result in interesting hysteretic effects. 
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1 Introduction 

The nonlinear diffusion equation : 


u = A (/(«)) (1) 

is the standard mathematical model for species diffusion in a continuous medium. The 
derivation of (jl| begins with the mass-balance equation: 

u + divh = 0, (2) 

which relates the time derivative of species concentration u to the divergence of the flux 
of diffusant h. Then, according to Fick’s law: 

h = -aVjtt, (3) 

the flux of diffusant is deemed proportional to the chemical-potential gradient V /j through 
a constant, positive mobility a. The combination of (j2| and ([3]) yields 

ii = aAg, (4) 
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and, for 'if'(-) the derivative of the coarse-grain free-energy mapping ^(-), the equation of 
state: 

p = i)' (u) (5) 

is enforced, so that ([TJ) is arrived at on setting: 

f(u) = aefi(u). (6) 

When modeling diffusion coupled with phase separation, one usually selects a non- 
convex coarse-grain free energy, the typical example being the double-well potential: 

if(u) = ku 2 (u — l) 2 , (7) 

with k a positive constant. In this case, ([TJ) is ill posed, and a regularization of the 
equation of state is in order, the most popular choice being the elliptic regularization-. 

H = — XAu, (8) 

with A a positive constant related to surface tension, leading to the celebrated Cahn- 
Hilliard equation pjj: 

u = aA — XAp). (9) 

An alternative approach was explored in [25], where ([8]) was replaced with 

p = + f3u, (10) 

with (3 a positive constant, whence the following viscous regularization of the diffusion 
equation: 

u — a/3Au = aA (11) 

which was shown in the same paper to be well posed when if'(-) is a locally-Lipshitz 
function fulfilling certain growth assumptions. 

Besides [25] . there is a substantial amount of mathematical literature devoted to the 
analysis of ( |TTj ) and of similar models. The case when is decreasing for large values 
of the argument has been considered in [33]. Related models have been studied in [2] and 
0 - 

The papers ra and [26] are concerned with the behavior of the solutions in the 
vanishing-viscosity limit /5 —^ 0 and provide insight on the hysteretic properties of the so- 
called entropy solutions. Such vanishing-viscosity approach provides a selection criterion 
for solutions outside the standard Sobolev setting in the degenerate case fd = 0, as shown 
in [2Ej and [32]. Forward-backward parabolic equations leading to hysteresis have also 
been considered in [35] . 

System obtained by combining both elliptic and viscous regularization have been 
considered in [T2] and [13], whose vanishing-viscosity limit has been studied in [31] . More 
sophisticated generalizations of the Cahn-Hilliard system that still incorporate a viscous 
constribution have been proposed and investigated in mmm- 

Despite the impressive amount of analytical literature concerning ([8]) , only a few 
references are available concerning its justification and interpretation. In [23 it was 
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shown that may be recovered as a suitable limit of the equations describing the 
motion of a mixture of two fluids. The possibility of giving a rational position to ([8]) is 
also intimated in [16] . A derivation based on microforce balance is also provided in [21] . 

Here, following the point of view of [18], we select a quite encompassing constitutive 
class (c/. (26)) and, within this class, we find the most general constitutive equations 
consistent with the second law of thermodynamics (c/. (36)), appropriate to the isother¬ 
mal context. In particular, we are able to retrieve (10) as a special case (c/. (37)). In 
addition we propose a proof of existence of weak solutions for a certain class of polynomial 
coarse-grain free energies. One of our motivations for presenting this proof is to illustrate 
how the viscous regularization affects the relevant estimates. 

The last section of this paper collects a handful of additional issues: we propose a non¬ 
smooth variant of (10) (c/. (89) and (90)), and we argue that this variant should produce 
interesting hysteresis effects; then, we cursorily examine singular free energies and we 
discuss the position of prescriptions like (10) within the non-standard thermodynamical 
setting put forth in [27]. 


2 Standard background 


The point of view promoted by Gurtin in H2, and illustrated in the recent monograph 
[18]. is that chemical potential must be treated as a primitive field representing the 
amount of energy carried by a unit amount of diffusant. In accordance with this notion, 
given any part V of the body where the diffusion process takes place, the quantity 



( 12 ) 


must be interpreted as the amount of chemical energy supplied to V per unit time. In the 
absence of mechanical interactions, the appropriate version of the dissipation inequality 
is: 


ipd x < T(fP). 


(13) 


V 


By combining (12) and (13), using the divergence theorem, and a localization argument, 
one obtains the inequality if + fi divh + h ■ V/i < 0; then with the help of the mass-balance 
equation (|TJ) , one easily arrives at the local dissipation inequality : 


ip — /Jit + h • V/i < 0. 


(14) 


It is standard practice in continuum mechanics to exploit dissipation inequalities to 

• single out those quantities that should be the object of constitutive specification: 
in the present case, the triplet (-0,/i, h); 

• determine what fields should appear as independent variables in these specifications: 
here, the triplet (u,u, g), with 

g = V/i. (15) 

• discard, in the manner of Coleman and Noll [7], “unphysical” constitutive choices, 
based on the requirement that the dissipation inequality be never violated along 
any realizable process. 
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An example of application of this procedure may be found in the recent monograph [T£] . 
whose Section 66 presents a fully-fledged constitutive theory for single-species transport 
coupled with elasticity. The starting point of the theory in question are certain provisional 
constitutive equations which, when strain is neglected, take the form: 

-0 = #u), = h = h(tt,g). (16) 


It is shown in [18j that consistency of (16) with the dissipation inequality (14) demands 
that the constitutive mappings delivering free energy and chemical potential be related 
by 


K u ) = v( u )^ 


and that 


h(u,g) = —M(u, g)g 

with the tensorial-mobility mapping M satisfying: 

M(u, g)g ■ g > 0. 


(17) 

(18) 

(19) 


It is important to notice that: 

• the relation ([5]) between chemical potential and free-energy mapping is recovered 
as a derived assertion; 


the (isotropic) Fick’s law 


can be recovered as a special case of (18), by choosing: 
M(tt,g)=al, (20) 


where the constant a is non-negative, as demanded by (19). 

For the sake of completeness, we briefly recapitulate the steps leading to @-@), and 
we refer to m for details. We begin by noticing that, on account of ( fl6| , the dissipation 
inequality (14) becomes: 


(w(u) - p(u))ii + h(tt, g) ■ g < 0. 


( 21 ) 


Following Coleman & Noll [7], we consider a process such that, at a given point, and at 
a given time, g = 0, with u and u attaining arbitrary values. At that particular point 
and time, the dissipation inequality specializes to: 


(ip'iu) — p{u))u < 0. 


( 22 ) 


We insists on asking that (22) be satisfied for whatever choice of u and u. Such require¬ 
ment can be met only if the constitutive mappings for free energy and chemical potential 
are related by (17). 

In view of (17), what is left with plj ) is the so-called residual inequality. 

h(w,g)-g<0. (23) 


Further conclusions can be drawn from (23) by fixing u and g, and by looking at the 
function: 

A i->- h(u, Ag) ■ g. (24) 
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By (23), the function specified in (24) changes its sign at A = 0. Moreover, this function 


is smooth, granted that the constitutive mapping h is smooth. Then, we have necessarily 
h(w, 0 ) ■ g = 0. Since u and g can be chosen arbitrarily, we conclude that: 

h(w, 0 ) = 0 . (25) 

In words: if the gradient of chemical potential is null, then the flux of diffusant is null 


as well. As a consequence of (25), the constitutive mapping h admits the representation 


(18), with (19) being required by the residual inequality (23). 


3 The parabolic regularization and its generaliza¬ 
tions 


Willing to explore constitutive dependencies more general than (16), we notice that the 


free-energy imbalance prompts the inclusion of ft in the set of independent variables 
appearing in the constitutive equations. We start from the following generalization of 


(16): 


0 = 0(w, ft), p = p(u, ft), h = h(it, ft, g). 


(26) 


We assume that the constitutive mappings appearing in (26) are all smooth. Then, the 


dissipation inequality (14) takes the form: 


(<9 u 0(w, ft, g) - n(u, ft))ft + h (u, ft, g) ■ g + du^(u, ft, g )u < 0. 


(27) 

If we assume that ( [27] ) holds true for whatever continuation of a process, and hence for 
any arbitrary assignment of it and g, we see that <9^0 = 0. Thus, there exists a function 
0 such that: 

0(u, ft) — 0(tt), (28) 


and (27) becomes: 


(<9 u 0(w) - n(u, ft))ft + h(tt, ft, g) ■ g < 0. (29) 

Arguing as in [31 Appendix], we set g = 0 in the above inequality, and we reckon that 
the function A (->• <9 u 0(w) — Jl(u, A) changes its sign at A = 0. As this function is smooth, 
we conclude that it vanishes for A = 0, and hence: 


K u , 0) = <9 u 0(w). 


It follows from (30) that exists f3(u,ii) such that 

/7(u, ft) = <9 M 0(w) + f3(u, u)u. 

Next, we define 


and 


h(«,g) = h(u,0,g), 

f dji(u, 0, g) if ft = 0, 

m(it, ft, g) = ^ h(w, ft, g) — h(w, g) 


otherwise. 


(30) 

(31) 

(32) 

(33) 


u 
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Then, we can write 


h(w, u, g) = h(u, g) + in (u, u, g )u. 


(34) 


On taking (31) and (34) into account, we arrive at the following form of the dissipation 
inequality: 

- f3(u, u)u + g • m(w, u, g)u + h(n, g) ■ g < 0. (35) 


On setting it — 0 and on insisting that 
the representation (18) for h. 


holds for whatever choice of g, we recover 


Now, on taking into account (28), (31), (18), and (34), we can rewrite the constitutive 
equations (26) as: 

4> = ip(u), n — d u tj}(u) + P(u,ii)ii, h = -M(m, g)g + iim(u,ii, g). (36) 


We are now in position to interpret the parabolic regularized equation (jll) as a conse- 

(37) 


quence of the following specifications in (36): 
P(u,u) = p, 


M(u, g) = al, m(tt, u, g) = 0, 


where the positivity of the constants a and P guarantees that the dissipation inequality 
is never violated. 


4 Existence of weak solutions in the viscous case 


In this section we consider the following system: 

it + div(aV/i) = 0, 
H — Pit + 


(38a) 


with a and P strictly positive constants. We obtain system (38a) by combining the mass 


balance equation Q with the constitutive prescription (j3j) and the viscous regularization 
(10). For fl C M 3 a domain with smooth boundary T, and for T > 0, we prove existence of 


weak solutions to (38a) in the parabolic domain Q — x (0, T) with the initial condition 


w(0) = u 0 in Q, 

and with the (possibly) non-homogeneous Neumann condition: 

Vp ■ n = h on S, 


(38b) 


(38c) 


imposed on the parabolic boundary E = T x (0, T), with n the outward unit normal on 
T. The treatment of the non-homogeneous Dirichlet condition for /i, can be found in [5]. 

The proof proposed by Novick-Cohen and Pego in [25j is based on regarding (38aJ) 
as an ODE in a suitable Banach space. Here we use a different approach, based on the 
Galerkin method and on compactness arguments (see for instance [19j or [29j). Our main 
purpose here is to illustrate where the viscous regularization comes in handy as far as 
existence of solutions is concerned. 
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For typographical convenience, we henceforth write 0 in place of 0; furthermore, we 
use the symbol Ci to denote a generic positive constant that depends on the index i: 
moreover, given p > 1 we denote by 


P = 


P 


p — 1 


(39) 


the conjugate Holder exponent of p; moreover, as a rule, we do not relabel subsequences. 
Following standard notation, for B a Banach space and r > 1 we denote by 1/(0, T; B) the 
Z/-Bochner space of H-valued functions defined on the interval (0, T), and by H 1 ( 0, T; B) 
the corresponding Sobolev-Bochner space. For the definition of these spaces, we make 
reference to Sections 1.5 and 7.1 of [25]. When elements of these functions spaces shall 
be involved in the definition of double integrals both with respect to the time variable t 
and space variable x, the dependence on the latter shall be left tacit. 

Assumptions. We assume that the free energy is twice continuously differentiable, and 
that its second derivative be bounded from below, namely: 


0"(r) > -M 0 


(40) 


for some Mq > 0, for all r G M. Moreover, we assume that there exist Mi > 0, i — 1... 5 
and p G [2,6) such that, for all r G M, 

- Mi + M 2 \r\ p < 0(r) < M 3 + M 4 |rf (41) 


and 

Finally, we assume: 


|0'(r)| < M 5 (l + \r\ v 1 ). 


h G 7/(0, T] L 2 (T)), 


(42) 

(43a) 


and 


u 0 G H\Q), 


(43b) 


with p the same exponent as in (41) and 
we are going prove the following: 


. Under the aforementioned assumptions, 


Theorem 1 (Existence of weak solutions to (38a)). There exist 


u G L°°(0, T; H\Q)) n H\ 0, T; L 2 (H)), 

peL 2 (o,T-H\n)), 


such that 



u(t)v + oV/x(f) • Vu) dx 
f3u + 0'(u) a.e. in Q 1 


u(0) = uq a.e. in 


h{t)v dr 


Vu G H 1 ( 0), for a.a. t G (0, T), 
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(44a) 

(44b) 

(44c) 





Proof. For the sake of readability, we split the proof into a sequence of intermediate steps. 


Step 1. Selection of a basis. We denote by {n ri }^ =1 the eigenfunctions of the Laplace 
operator with Neumann boundary conditions: 

—Av n = X n v n in fl, (45a) 

Vv n • n = 0 on T, (45b) 


and we define 


V n = span(ui, ...v n ). 


(46) 


The set U n V n is dense in // 1 (h2) and in L P (Q), where p is the exponent in (41) and 
(42). Without any loss of generality, we can assume that the collection {r; n }^ =1 is an 
orthonormal system for L 2 (Q), that is: 


XnVk d.C dnfc; 


(47) 


m 


and that the first element of the basis is the constant function: fi(a;) = 1/yfOj. We 
approximate the initial datum u 0 through a sequence {wo.nl^Li of functions having the 
form 


U 0 , n (x) = ^2 


^0 .n.i'ViyX) 


(48) 


i=l 


such that 


Uq < n u o strongly in d/ 1 (h2), (49) 

with convergence holding also in L P (Q) because of the Sobolev embedding. Note that, 


thanks to our assumption (43b) on the initial datum and to our growth assumption (41) 


on the coarse-grain free energy 'if, we have 

if{u o, n ) if(u 0 ) strongly in 11(0). 
Thus, in particular, there exists a positive constant Cq such that: 


(50) 


V : (ui, v )(U < Co- 


(51) 


Step 2. Construction of an approximating sequence. We introduce the Galerkin ap¬ 
proximations 


n 

U n (x,t) = ^2a nti (t)vi(x), (52a) 

2—1 
n 

Unix, t) = ^2 bn,i(t)vi(x), (52b) 

2=1 
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and we look for time-dependent coefficients (a n}i (t), b nji (t)) i= i___ n such that 
/ (u n v + aV fi n ■ Vn) dx — I hv dr for all v G V n , 


(53a) 
(53b) 

(54) 

To this aim, by testing (53a) and (53b) by v t , we find that the system (53) is equivalent 
to: 


at all times, and 


fi n v dx = / {/3u n + ^{u n ))vdx = 0 for all v € V n , 

\ Jn 

U n ( d ) Wo,n- 


®n,i(f) T a\ib n ,i(t) -/Tj(f), 


i = 1, • • • ,n, 


(55) 


bn,lit') /3hri,i(^) T G n j, (®n,l (f) j • • • Q j n,n(t')') i 

with initial conditions: 

®n,i(0) ®0 ,n,ii * — 1, • • • , 77>, (56) 

where ifj(i) = f r h(x, t)vi(x)dr and where the functions : M n —» M are defined 

by G n>i (ri,... , r n ) = r j v j( x )) v i ( x ) 6.x'. Since the eigenvalues A, are non 

negative, the two groups of equations in ( |55j ) can be combined to obtain a n x n system 
of ordinary differential equations: 

(xXiGn^i (dr) ,l (t ), • • • fln,n(f)) 


b>n,i {t ) 


1 + a/3Xi 


i = 1... n, 


(57) 


which together with the initial conditions (56) defines a Cauchy problem. 

By the smoothness of ^ and of the basis functions, the functions G n ^ are smooth as 
well, and for each i e 1 ...n the right-hand side in (57) is locally Lipschitz continuous. 
Thus, according to standard ODE theory, the Cauchy problem (56)—(57) has a unique 
solution on a non-empty interval (0, T n ). If T n < T , we must show that the solution admits 
a continuation up to time T. We achieve this goal in the next step by deriving bounds on 
the coefficients in (52) based on an estimate which mimicks the natural energetic estimate 
for the original system (38a). 

Step 3. Energetic estimate. We consider a subinterval (0, t) of (0 ,T n ), and for each 
s G (0, t) we test (53a) and (53b) with p n (s) and —u n (s), respectively. On adding the 
resulting equations and on integrating over (0, t) we obtain: 


ip(u n (t))dx+ / / (a\V/j n (s)\ 2 +/3ul(s)) dxds 


o Jn 


h(s)n n (s) dTds + / , 0( M o,n)da:. (58) 


o Jr 


Thanks to the coercivity assumption in (41), and thanks to (51), we have 

ft 


IQ 


M 2 / \u(t)\ p dx + / / (a\Vfj n (s)\ 2 + pu 2 n (s))dxds 


o Jn 


< Mi|f2| + Co 

cT 



h(s)p n (s)dTds, (59) 


o Jr 
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where M\ and U 0 are the constants in, respectively, (41) and (51). Using, in the order, 


the trace theorem, Poincare’s inequality, and Young’s inequality, we have, at each time 
(whose specification we omit to keep our notation terse), 


I hn n dr < 11 11L 2 (r) 11 l^n 11 L 2 (r) 

< Cl\\h\\ L ^{Y)\\^n\\ H i{Q.) 

< C 3 \\h\\ L 2 {r) ( ||V/X n ||z, 2 (n) + / Pndx 


in 


<C 4 (-||hr i2(r) +^(||V^rn' 2( o) + 


/in d.X 


(60) 


where we observe that (recall (39) and the assumption on p in the sentence preceding 


©): 

p' < 2. (61) 

On choosing v = 1 as test in (53b), we obtain f Q p n dx = f n ( aii n + i/>'(u n ))dx. Then, on 
taking the growth assumption (42) into account, we get f n p n dx < C^l + J Q (|tt n + 

|u n | p_1 ) dxj, and hence: 


/in dx 


in 


<C 6 [ 1+ / (\u n \ p ' + \u n \ p ) dx 
Jn 


(62) 


Thus, (60) and (62), along with our assumption (43a) on h, and (61), yield, for every 
5 > 0, 


hp n dT < ^\\h\\ p L2(r) +S p 'C 7 (l + ||V/i n ||J a(n) + / (|M„K + K| p )da;) 


<^ + 5 p 'C 9 ( 1+ / (\Vp n \ 2 + \u n \ 2 + \u n \P)dxy 


(63) 


By gathering (59) and (63), we obtain 


M 2 / \u n {t)\ p dx + / / (a\Vp n (s)\ 2 + fJu 2 n (s)) dxds 
Jn JoJn 

< Ci + § + S P 'C 9 [ f (|V/i n (s)| 2 + |u„(s)| 2 + |u n (s)| p )da;ds, (64) 
Jo Jn 


Sp 


whence, on taking S sufficiently small, 

M 2 [ \u n (t)\ p dx + ^ f f (a\Vp n (s)\ 2 + (3iil(s))dxds < C w (l + ( f \u n (s)\ p dxds). 
Jn ^ Jo Jn v JoJn ' 

( 65 ) 

By using Gronwall’s lemma in (65), we obtain that the coefficients a[ n \t) are bounded 


for t G (0, T n ) by a constant that does not depend on n. It then follows that the solution 


of (57) can be continued up to time T and that the bounds: 

|| w n|U°°(0,T;LP(fi)) < Cll, 

||||L 2 (Q) < c 12 , 


(66a) 

(66b) 
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hold uniformly with respect to n. Moreover, we have that || V/x n ||x, 2 (Q) is uniformly 
bounded, a fact that, together with (62) and (66) yields 


11 h"n 11 L 2 (OfT-,!! 1 (f2)) < Ci3- 


(67) 


Step 4- Limit passage in the diffusion equation. The estimates (66b) and (67) are 


what we need to pass to the limit in the linear equation (53a) to obtain (44a). Indeed, 


thanks to these estimates there exists a subsequence (not relabeled) such that 

u n —y u weakly in i7 1 (0,T; L 2 (0)), 
p, n —* fi weakly in L 2 (0, T; H 1 (12)). 


( 68 ) 

(69) 


Now, given any k < n, we have that Vk is a legal test in ( |53a ). Thus, for every for every 
t G (0, T ) we can consider: 


° = s 




(u n (s)vk + aVfin(s) ■ Vvk)dxds = 0 Vn > k, 


(70) 


for all 0 < h < T — t. On letting first n —> oo and then h —> 0 in (pOk the Steklov average 


in (70) converges for a.e. t G (0,T): 


u{t)v k + aS7pi{t) ■ Wv k = 0 Vt G (0, T) \ E k , 


(71) 


where meas (Ek) = 0 is a set that depends on the particular k. Now, given that (72) holds 
for k arbitraty, we conclude that 


( u(t)v + aV n(t) ■ Vu)da; = 0 Vr: G UfeI4, Vf G (0,T) \ U k E k . 


(72) 


Since U^I4 is dense in 7/ 1 (0), and since the set U k E k has zero measure, we conclude that 


(44a) holds true. 


Step 5. Gradient estimate on concentration. In order to derive (44b) by passing to the 


limit in (53b) weak convergence does not suffice because of the possibly nonlinear term 


f>'(u n ). We therefore seek extra compactness properties for u n: a piece of information 
that we shall gather from an additional estimate. In order to get this estimate, we first 


choose u n as tests in (53a) to get, for every s G [0, T], 


7 n (s)u n (s) + aVn n (s) ■ Vu n (s)) dx= h(s)u n (s) dx; 


(73a) 


then, we take v = —aAu n in (53b) and we resort to by-parts integration to arrive at: 


a / V/x„(s) • Vu n (s)dx = a / (/3Vm„(s) • V«„(s) + (u n (s))\ Vw n (s)| 2 )dx, (73b) 

Jn Jn 
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which again holds for all s G [0, T]; next, we subtract (73b) from (73a) and we integrate 
with respect to s on (0, t) to get: 


2 / (ul(t) + a/3\Vu n (t)\ 2 )dx + a \ / if"(u n (s))\ Vw n (s)| 2 dxds 


rt 

'o Jn 


1 

2 


(u 2 {x, 0) + a/3\X7u n (x, 0)| 2 )dx + / / h(s)u n (s) dTds (74) 

Jo Jr 


It is precisely at this point that the viscous regularization comes in handy. Indeed, with 
the aid of Assumption (40) we obtain from (74) the following inequality: 


1 

2 




u r 


(t) + a/3\Wu n (t)\ 2 )dx < aM 0 / / | Vu n (s)| 2 ds + / / h(s)u n (s) dTds 


o Jn 


o Jr 


1 

2 


u 2 n {x, 0) + a/3\ V« n (x, 0)| 2 )dx. (75) 


Next, by making use of (43a) and by invoking the trace theorem, we estimate: 


h(s)u n (s) drds < - 


/r (s) drds + - 


w;(s) drds 


o Jr 


< C i4 (1 + / ||w n (s)||tfi(n) ds ] • 


'0 


(76) 


On combining the inequalities (77) and (76), and by recalling (49), we arrive at: 

rt 


\ u n{t)\\ 2 m(n) < C is ( 1 + 


K( s )ll/fi ( n) ds 


Now, the application of Gronwall’s lemma to (77) yields: 

< Cie, 


(77) 


(78) 


a bound that we are going to use in the next step to pass to the limit in the nonlinear 
equation (53b). 

Step 6. Limit passage in the equation for chemical potential. Since p < 6 and O is 
a three-dimensional domain, i/ 1 (0) is compactly embedded in L p (Ll). Consequently, we 
can use the bounds (78) and (66b), along with the Aubin-Lions lemma in the form of 
Corollary 4 of [30J, to establish the following convergence result: 

u n —y u in C([0,T];L p (0)). (79) 


We note on passing that (79), along with with (49), entails the initial condition (44c). 

Now, on rewriting Assumption (42) as \if'(r)\ < Af 5 (1 + |r| p//p ), we immediately see 
that standard results concerning mappings between Lebesgue spaces (see for instance 
Theorem 1.27 in [29]) do apply, and we can conclude that the Nemytskii mapping u i —> 
fj’(u) is bounded and strongly continuous from L p {Ll) to L p (Q). Thus, the convergence 


in C([0, T]; L p '(Q)) 


(80) 
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follows from (79). Now, pick arbitrary functions a G C([0, T]) and <p £ E p (0), and let 
cp n G Vn be a sequence of linear combinations of basis functions such that 


(p n —y (p strongly in L p {fiV). 


(81) 


On testing (53b) by a(t)cp n at each particular time level t G (0, T) and on integrating 
over (0, T ) we get: 



o Jn 


~ 0u n (t) + ip'(u n (t)))a(t)<t> n Axdt = 0. 


(82) 


By relying on the weak convergences (68) and (69), and on the strong convergences (80) 


and (81), we can pass to the limit in ( [82] ) to get: 

r T 


(p(t) — fiu{t ) — ip' (u(t))^) cp dx) a(t)dt = 0 


(83) 


Jo v Jn 

By the arbitrariness of the function t t —> a(t) we conclude that 


in 


(/r(f) — /3u(t) — ip'(u(t))^(p({x for a.e. t G (0, T). 


(84) 


Eventually, the arbitrariness of <p in (84) entails (44b). 


Step 7. Verification of the initial condition. In order to show that the initial condition 


(44c) is satisfied, we begin by observing that, at the approximation level n, 

n ii n (pdxdt = — / / u n (pdxdt— / u o ,n0(O)da; 

i Jo Jn Jn 


(85) 


for every smooth test function <p G C°°(Q) such that <p{T) = 0. By passing to the limit 


in (85), using the weak convergence of the sequences in (69) and the convergence of the 


initial conditions for the stated in (49), we arrive at: 



ucpdxdt = — / / ucpdxdt— / u o 0(O)da:. (86) 

i Jo Jn Jn 

On the other hand, since u G L 2 (Q ), we obtain, through integration by parts, 

n u(pdxdt = — f ( ucpdxdt— f w(O)0(O)da:, (87) 

i Jo Jn Jn 


whence, by a comparison with (86), 


/ u(O)0(O)dx = / u 0 cp( 0)dx, 

In Jn 


and thence the initial condition (44c), by the arbitrariness of </>(0) 


( 88 ) 

□ 
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5 Additional remarks 


Non—smooth dissipation potentials and hysteresis. Our discussion in Section 3 
suggests that the class of viscous regularizations described by ( 10 ) may be replaced with 
more encompassing prescriptions for /i, provided that the result is consistent with the 
dissipation inequality. For example, one may want to drop the smoothness assumption 
we have made so far and consider the inclusion: 


A* - V («) 6 9?(“) 


(89) 


with C a proper, lower-semicontinuous, convex dissipation potential. We here restrict 
attention to the special case: 

CM = ^r 2 + 7 |r|, (90) 

with f3 > 0 a viscous-damping parameter and 7 > 0 a threshold parameter. In this case, 
the subdifferential set of ( at r is given by: 


d((r) = 



[-7, + 7 ] 


if r ^ 0 , 
if r = 0 . 


(91) 


Accodingly, the inclusion (89) is equivalent to: 


u 


fi = 0 '( u ) + (3u + 7 -pry if ii 7 ^ 0 , 


m 


v-i' (u) e [-7,7] 


if u = 0 . 


(92a) 

(92b) 


The prescriptions (92), when combined with the PDE that results from (J2| and ([ 3 ]) , 
namely, 

u = aA/i, (93) 

together with appropriate boundary conditions, yield an evolution problem in the un¬ 
knowns u and p. The analysis of this problem, which is more involved than the case of 
viscous regularization, shall be presented in pj. 

We expect that the solutions of (92) (93) display hysteretic behavior , even if the free 
energy mapping -0 is convex. As an example to support this point, we consider that the 
free-energy mapping is quadratic: 


i>(r) = jAr 2 , 


(94) 


so that (92) becomes: 


u 

Ittl 


p - Ku G [- 7 , 7 ] 

We impose a homogeneous initial condition: 

u(x, 0 ) = 0 f( 


if ii 7 ^ 0 , 

(95a) 

if u = 0 . 

(95b) 

x G fh 

(96) 


14 








Moreover, instead of prescribing the flux at the boundary as in Section 4, we impose the 
Dirichlet condition: 

p(x,t) = Az(^-^j for all x G T, (97) 

where A > 0 is the half-amplitude, r > 0 is the period and z(s) the zigzag function ( cf. 
Fig l.(a) below) defined on [0,1) by: 

( 4s if s G [0, 1/4), 

z(s) = < 2 -4s if s G [1/4, 3/4), (98) 

{ -4 +4s if s G [3/4,1), 


and extended periodically to M. 

The particular form of the boundary datum prompts us to rewrite (92) and (93) using 
the dimensionless time : 

s — — (99) 

r 

in place of t as independent variable. Accordingly, all fields are henceforth functions of x 
and s, and a superimposed dot denotes the derivative with respect to s. Then, (92) and 
(93) become, respectively, 


„ (3 . u 

H — Ku-\ - U + 7—r 


\u\ 


h - Ku G [- 7 , 7 ] 


if it 7^ 0, 
if u = 0, 


and 


u 


— = aA/i. 


r 


We are interested in the (formal) limit 


t —y 00 , 


(100a) 

(100b) 

( 101 ) 

( 102 ) 


which describes a regime when the boundary datum changes slowly. I11 this regime, ( 100[ ) 
and (101) become, respectively, 


and 


„ u 
p, — Ku + 7,—r 
\u\ 

p-Ku G [- 7 , 7 ] 

aAp = 0. 


if it 7^ 0, 

(103a) 

if u = 0, 

(103b) 


(104) 


Now, in view of the boundary condition (97), it follows from (104) that // is spatially 
constant in at each time, that is: 


p(x, s) = Az(s). (105) 

Given that the initial condition for u does not depend on x as well, we can look for 
constant-in-space solutions-. 

u(x,s) — v(s), (106) 


15 












with i?(s) satisfying 


Az(s ) = Kv(s ) + 7 
Az(s) — Kv(s ) G [- 


1 ; 0 ) 

|h(s)| 

■7,7] 


ifh( S )^0, 
if v(s) = 0 . 


(107a) 

(107b) 


Given that z(0) = 0 by (98), system ( |107 ) has a unique solution v G H 1 ( 0, T; M) satisfying 
the initial condition u(0) = 0 |3U Lemma 2.9]. This solution, which in general is described 
by a play operator £ such that v(-) = £(Az(-), r>(0)), can be worked out explicitly in the 
present case. 

If the half-amplitude A of the oscillation is smaller than or equal to the threshold 7 , 
the solution is null at all times; otherwise, for V the projection on the closed interval 
[—A + 7 , A — 7 ], the solution is given by: 


Kv{s ) = 


0 if s G [0, 2a) > 

V(Az(s-2a)) if s G [£4, + 00 ) , 

as shown in Fig. l.(a) below. The graphs of Az(s) along with Kv(s ) are shown in Fig. 



Figure 1: (a) plots of Az(s) and Kv(s ); (b) parametric plot of s (Az(s), Kv(s)). Here 
s — t/r is the rescaled time. 


2. In particular, hysteresis is apparent from Fig. l.(b). 


Singular free energies. One may want to consider the following singular free energy 
suggested by the regular solution model [1]. In this model, the entropic contribution 

2> e {u) = kr logr, k > 0 , (108) 

to the free energy is augmented by a quadratic enthalpic term through an interaction 
parameter x > 0 , so that the total free energy is: 


2>{r) = V’e( r ) + "0e(l — r) + yr(l — r) if r G (0,1), 


(109) 


with 2’{ r ) = +°° otherwise. This prototypical energy is, for instance, relevant when 
describing solvent diffusion in polymer gels (see, e.g., manaro. 

Now, the issue with the singular terms in (109) is that, at variance with energies having 
polynomial growth, ^'(r) is not controlled by V’e( r )- In fact, as r G (0,1) approaches 0 
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Figure 2: Plots of the free energy in (109) and of its derivative (respectively, solid and 
dashed line). 


or 1, we have that Vy(r) tends to null, whereas |^'(r)| blows up, as shown in Fig. 2 


below. In this case, the argument leading from (58) to the energetic estimates (66) 


breaks down, even when carried out formally. Indeed, the bound (62), which is essential 


to that argument, follows from (the discrete counterpart of) the equation 

p = f3u + 


( 110 ) 


whose exploitation is contingent on having a control on For energies with polyno¬ 
mial growth, such control is guaranteed from i/)'(u ) being bounded by but this is 


not the case when ip(u) has the form (109). 

The fact that a basic energetic estimate fails in this situation is the symptom of a 
structural feature that would be fatal to any naive attempt to handle the problem by 
approximation of the free energy. Suppose in fact that, by mimicking m , we approximate 
-0 e with the following regularized functions: 


if. ; e (r) if r > e, 

lM r ) = { k /r 2 

kr\oge+-y - e) if r<e, 


( 111 ) 


and we consider a e-parametrized sequence of problems with the following regularized 
free energy: 

A(r) = ijj e ,e(r) + ^e,e(l ~ r) + xr(l - r). (112) 


The energies specihed by (112) are smooth and have quadratic growth as |r| —>• +oo. 


Thus, Theorem [l] can be applied to establish the existence, for each particular e, of an 
approximation 


u, € L°°(0,T-,H\n)) n// 1 (0,T;L 2 (SJ)), 


(113) 

(114) 


that satisfies 


/ u £ {t)v + aS7fi £ (t) ■ Vv = / hv Vu e Ff 1 (h2) for a.a. t 6 (0,T), 

Jn Jr 

fi £ = au £ + ip' E (u £ ) a.e. in Q , 
w e (0) = uq a.e. in Q. 
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(115a) 

(115b) 

(115c) 










Since approximates f) in the sense that 


V4(r) —>■ f>'(r) as £ — y 0 for all r e (0,1), 


(116) 


one may hope that, as £ —> 0, the solution of (115) converge to a limit, and this limit 
satisfy the original problem. 

In order to establish such convergence, however, one should rely on estimates analo¬ 


gous to (66), (67), and (78). The problem however is that the bounds in these estimates 


would not be uniform as e —> 0, because if £ (r) bounds -0'(r) on ly through a multiplicative 
constant that blows up as e —» 0. 


On the other hand, if the Neumann boundary condition (38c) were to be replaced by 
the Dirichlet condition: 

H — H b on T, (117) 


as done in [5], then an energy balance similar to (58) would hold. In this case, it would be 


possible to recover the estimates (66) on /i and u without recourse to (62), using instead 


the control on /i provided by its trace. Eventually, one would be able to control ip'(u) by 


comparison in (110). 


All in all, if the boundary condition (117) is enforced, the estimates 


IMU oo (0,T;.H' 1 (fi) < C, 
l“ £ || l 2 {q) < C, 


(118) 

(119) 

( 120 ) 


11 /A 11L 2 (0,T 1 (f2)) < C, 

hold uniformly with respect to e. Using these estimates, we can pass to the limit in 


(115b) by replicating the argument in [Hi- 

Alternative frameworks. The point of view put forth in |TS3, which we illustrate 
in Section 2, appears to be less general than [16], where a microforce balance: 


div£ + 7T + 7 = 0 (121) 

augments the mass balance (J2| and chemical potential is included in the list of inde¬ 
pendent state variables. In particular, the internal microforce w, rather than chemical 
potential, is the object of a constitutive prescription (see also the discussion in [5J). How¬ 
ever, since we do not consider surface-tension effects, in our case the vectorial microstress 
vanishes: 

£ = 0 , ( 122 ) 

and, given that the external microforce 7 vanishes in any actual evolution process: 


7 = 0, (123) 

the microscopic force balance (121) reduces to the statement tt = 0, a statement that, 


when combined with the constitutive equation 

7T = //-'0 , ( W ), 


( 124 ) 


yields the equation of state ([5]). 
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An alternative framework for processes of phase diffusion, which hinges on the analogy 
between chemical potential and coldness as measures of orderliness , has been proposed in 


In that framework, the microforce balance (121) is retained, but balance of micro¬ 


scopic energy and imbalance of microscopic entropy replace, respectively, mass balance 
and dissipation inequality. Consistent with this point of view, the microscopic entropy 
flux, which is involved in the microscopic entropy imbalance, is assumed to be propor¬ 
tional through /i to the microscopic energy flux, which appears in the microscopic energy 
balance. With this assumption, the microscopic entropy imbalance can be combined with 
the microscopic balances of force and energy to arrive at the following counterpart of the 
conventional reduced dissipation inequality : 


if < )' + p, n • V/i — 7rii + £ • Vii, 


(125) 


which involves, besides the aforementioned flux of microscopic energy h, a microscopic 
free energy if and a microscopic entropy r). In particular a collection of constitutive 
prescriptions are sorted out, leading to diffusion models a la Cahn-Hilliard which have 
been studied in 0E]. These prescriptions include: 

if — Vu,h), 7 t = -d u if{u,Vu,n), £ = d^ u fl>(u, Vu, /i). (126) 

In particular, if the constitutive mapping delivering the free energy density has the form: 


if(u, Vu,/i) = — pu + if(u) + -ft|Vu| 2 , 


(127) 


then the microforce balance (121) with 7 = 0 leads to the equation of state 


In this framework, the analogue of the viscous regularization (10) would be: 


Moreover, the inclusion 


7 T = —d u fll(u, Vu, fl) — flu. 


- IT - 9„V>(m, Vu, n) 6 9C(*)> 


(128) 


(129) 


would take the place (89). The analytical consequences of both (128) and (129) are still 
to be explored. 
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